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Plane partitions naturally appear in many problems of statistical physics and 
quantum field theory, for instance, in the theory of faceted crystals and of topolog- 
ical strings on Calabi-Yau threefolds. In this paper a connection is made between 
the exactly solvable model with the boson dynamical variables and a problem of 
enumeration of boxed plane partitions - three dimensional Young diagrams placed 
into a box of a finite size. The correlation functions of the boson model may be con- 
sidered as the generating functionals of the Young diagrams with the fixed heights 
of its certain columns. The evaluation of the correlation functions is based on the 
Yang-Baxter algebra. The analytical answers are obtained in terms of determinants 
and they can also be expressed through the Schur functions. 

I. INTRODUCTION 

The theory of plane partitions is a classical chapter in combinatorics (lL Statistics of 
plane partitions with respect to natural probabilistic measures was studied in 2^ Q]. Locally, 
plane partitions are equivalent to a tiling of a plane by rhombi for lozenges). Among many 
results obtained in this direction we should mention the paper |4[ where the phenomenon of 
the "arctic circle" was described for boxed plane partitions. In the paper |5| the variational 
principle applicable to a variety of such problems was developed. Correlation functions for 
random plane partitions were studied in For boxed plane partitions distributed 

uniformly they were computed in the bulk of the limit shape (inside of the arctic circle) in 
for unrestricted plane partitions distributed as g' 71 "', where \tt\ is a volume of 



, and in 
a partition. 



Plane partitions naturally appear in many problems of statistical physics, for instance, in 
the theory of faceted crystals 7( and of direct percolation |s|] . Quite recently it was argued 
that there was a connection between topological strings on Calabi-Yau threefolds and crystal 
melting 
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In this paper we demonstrate the connection of a certain integrable boson type model 
and boxed plane partitions - the plane partitions placed into a box of a finite size. The 
natural dynamical variables in this boson model are the q— bosons. The algebra of qr-bosons 

n 

jlllj appear naturally within the quantum algebra formalism [10J. 

Boxed plane partitions are related to the special case of q— bosons, namely to a limit when 
the deformation parameter q tends to zero what corresponds to an infinite value of coupling 
constant of the q— boson model. In this special limit q— bosons are known as exponential 
phase operators of quantum non-linear optics. Notice that this is also the famous crystal 
limit for quantum groups jl^ . [l^j ]. 

Our analysis is based on the Quantum Inverse Scattering Method (QISM) and on the 
algebraic approach to the calculation of the correlation functions developed within this 



method ll 
mo 



15 



, . We shall show that the scalar product of the state vectors of the phase 

rrnm 

del [16( , |17[ , [18|] is related to MacMahon enumeration formula for boxed plane partitions. 



There is a well established connection between the theory of plane partitions and the theory 
of random processes. From that point of view the scalar product of a phase model is a 
generator of a point fermion-like random field jjj. The systematic application of the QISM 
allows to calculate different correlation functions appearing in the theory of the boxed plane 
partitions. 

II. g-BOSONS 

The q— boson algebra is defined by three independent operators B, and N satisfying 
commutation relations 

[B,B*\ = q 2N , 

[AT, St] = B\ [N,B] = -B, (1) 

and q, a c- number, is taken to be q = e -7 . We consider real 7 > 0. The q— boson algebra (JIJ 
has the representation in the Fock space formed from the g-boson normalized states \n) 

B ] \n) = [n + l]3|n+ 1) , B\n) = [n]*|n - 1), (2) 

where the "box" is 

1 - a 2n 



The integer numbers n > are called occupation numbers or the number of particles in a 
state \n): 

N\n) = n\n). (4) 
If q = 1 (7 = 0), the q— bosons become ordinary bosons, 

S->6,S t ->6 t ,iV=6 t 6; [6, 6 f ] = 1. (5) 

In the limit q — > (7 — > 00) the operators 5, i?t transform into the operators 0, 0' defined 
by the commutation relations 

[7V,0] = -0, [iV,0t] = 0t 5 [0,0t ] = 7r (6) 

in which it is the vacuum projector n = |0)(0|. The Fock states \n) can be created from the 
vacuum state |0) by operating by the phase operators, 

\n) = (<^) n |0), N\n) = n\n), (7) 

and 

^\n) = \n + l), </>\n) = \n-l), <f>\0) = 0. (8) 

One may verify that (j) an d $ can be expressed in terms of the Fock states, \n), as 

00 00 
(j) = \n)(n + 1|, (ft = ^2 \ n + 

n=0 n=0 

The introduced operator is "one-sided unitary" or an isometric, although 

= 1, 

one has 



= 1 - 

The operators © may be expressed in terms of ordinary bosons (JSJ): 

(j)= (tfb+iyh, $ = 6 + (6 t 6+ 1) _ 3. 



III. INTEGRABLE PHASE MODEL 



The phase model is a special limit of the integrable g-boson model 
by the L— operator [3], 



is defined 



/ 



LJu) = 



a n (u) b n (u) 
\ c n {u) d n {u) J y<p n u 

where parameter u G C, and cf) n , (j)^ are the operators (jBJ) satisfying commutation relations 



(9) 



On the local Fock vectors 

<f>j\0)j = o, 

The operator valued matrix (jHJ) satisfies the intertwining relation 

#0, v) {L n (u) <g> L n (») = {L n (v) <g> L n (u)) R(u, v) 

in which R(u, v) is the 4x4 matrix with the non-zero elements equal to 

R u (u,v) = R u (u,v) = f{v,u), 
R 22 {u,v) = R 33 (u,v) = g{v,u), 



(10) 



(11) 



(12) 



(13) 



R2?,(U,V) 



with 



f(v,u) 



u . . uv 

g{v,u) 



u 2 — V 2 



u 2 — V 2 



(14) 



Symbol (g> denotes the tensor product of matrices: (A <8> B)^ kl = AijB u . The R- matrix (fTTfl) 
satisfies the Yang-Baxter equation 



(7 <g> R(u, v)) (R(u, w) ® I) (I ® R(v, w)) 
= (R(v, w) ® J) (I ® R(u, w)) (R(u, v) ® I) 

The monodromy matrix is introduced as 

' A(u) B(u) 
C(u) D(u) 



T(u) = L M {u)L M ^{u)...L G {u) 



(15) 



(16) 



V 
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Matrix elements of this matrix act in the Fock space spanned on the state vectors 

M M 

i™> = ni"i>;=n(0]Ho>, z^=n, (I?) 

3=0 3=0 

where 

M 

io> = n io>i (is) 

i=o 

is the vacuum vector (generating state), and k( n m\ n i)j — ^im^kj- 

The commutation relations of the matrix elements of the monodromy matrix are given 
by the R— matrix (fT^j) 

R(u, v) (T(u) ® T(v)) = (T(v) ® T(w)) u). (19) 

The most important relations are 

C(u)B(v) = g(u,v) {A(u)D(v) - A(v)D(u)} , (20) 

C{u)A{v) = f(v,u)A(v)C(u)+g(u,v)A(u)C(v), (21) 

D(u)B(v) = f{v,u)B(v)D(u)+g(u,v)B(u)D{v), (22) 

[B(u),B(v)] = [C(u),C(v)]=0. (23) 

The relation (|19|1 means that the transfer matrix t(u) = trT(u) = A{u) + D(u) is the 
generating function of the integrals of motion: [t(u), t(v)] = for all u,v G C. 
The L— operator JSJ) satisfies the relation 

e< N "L n (u)ei<"' =ei<"'L„(u)e< N -,a'=\ 1 ° , (24) 

where C G C, and iV n is the number operator pip. From this equation and the definition of 
the monodromy matrix (fH)|) it follows that 

e^T(u)e^ z = e^ z T(u)e^, (25) 

where 

M 

N=J2 N n (26) 

n=0 



is a total number operator, and a z is the Pauli matrix (|24j) . The equation (|25p is equivalent 
to 

iVB(«) = B(u) (n + 1), (27) 
iVC(n) = C(u) (n-i). 
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(b) (c) (a) (d) 



FIG. 1: Vertex representation of the matrix elements of the L-operator. 

It means that the operator B(u) is a creation operator, while C(u) is annihilation one. 

The generating vector |0) (|18|l is annihilated by C(u) operator and is the eigenvector of 
A(u) and D(u) : 

C(u)\0) = 0, (28) 
A(u)\0) = a(u)\0), D(u)\0) = d(u)\0), 

with the eigenvalues a(u) = u^^ M+1 ^ and d(u) = u M+1 respectively. The iV— particle vectors 
are taken to be of the form 

N 

\q N ( Ul ,u 2 ,...,u N )) = Y[B( Uj )\0), (29) 

3=1 

and 

N\^ N (u 1 ,u 2 , ...,«jv)) = N\^ N (ui,u 2 , ..,%)). 

The vacuum state ()18)) is similar to the highest-weight vector in the theory of representations 
of Lie algebras. The state conjugated to (jSHj) is 

N 

(y N { Ul ,u 2 ,...,u N )\ = (OlUCiuj). (30) 

The dual vacuum (0| = |0) f , (0|0) = 1. It is easy to verify that (0\B(u) = and (0\A(u) = 
a{u){0\, (0\D(u) = d(u)(0\. 

One can visualize matrix elements of the L— operator as a vertex with the attached 
arrows (see FIG. 1). The matrix element b n (u) = ^ corresponds to a vertex (b), c n (u) = <p n 
corresponds to a vertex (c), a n (u) = u^ 1 to a vertex (a) and d n {u) = u to a vertex (d) 
respectively. 

Matrix elements of the monodromy matrix (JTfij) are expressed then as sums over all possi- 
ble configurations of arrows with different boundary conditions on a one-dimensional lattice 
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with M + 1 sites (see FIG. 2). Namely, operator B(u) corresponds to the boundary condi- 
tions when arrows on the top and bottom of the lattice are pointing outward (configuration 
(B)). Operator C(u) corresponds to the boundary conditions when arrows on the top and 
bottom of the lattice are pointing outward (configuration (C)). Operators A(u) and D(u) 
correspond to the boundary conditions when arrows on the top and bottom of the lattice 
are pointing up and down respectively (configurations (A) and (D)): 



i=0 



(B) (C) (A) (D) 



FIG. 2: Matrix elements of the monodromy matrix T(u). 



For example an operator B(u) on the lattice of three sites (M = 2) is: B(u) 
6q0!02 + u 2 4>\ + 4>\, and may be represented in a form (see FIG. 3). 



^u- 2 + 



B(u)-- 



+ i 



i=2 
i=l 
i=0 



FIG. 3: Vertex representation of the operator B(u). 



IV. SCALAR PRODUCTS AND PLANE PARTITIONS 

Recall that a partition A = (Ai, A 2 , ■■■,Xn) is a non-increasing sequence of non-negative 
integers Ai > A 2 > ... > A^ called the parts of A. The sum of the parts of A is denoted by 
|A|. 

A plane partition is an array ti^ of non-negative integers that are non-increasing as 
functions of both i and j = 1,2, ...). The integers it^ are called the parts of the plane 
partition, and | vr | = J^^ij is its volume. The plane partitions are often interpreted as stacks 
of cubes (three-dimensional Young diagrams). The height of a stack with coordinates ij is 
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7Tjj. If we have % < r, j < s and rc^ < t for all cubes of the plane partition, it is said that the 
plane partition is contained in a box with side lengths r, s, t. The symmetric plane partition 
is the plane partition for which 7Ty = itji. In (FIG. 4) the diagram corresponding to the 
boxed plane partition 4x4x5 



7T 



^5 3 3 2 ^ 

5 2 11 

4 110 

3 110 



(31) 



is represented. 




FIG. 4: A plane partition with a gradient lines. 

A plane partition inarxsxt box is equivalent to a lozenge tiling of an (r, s, t)-semiregular 
hexagon. The term lozenge refers to a unit rhombi with angles of | and To each plane 




a c d 

FIG. 5: Three types of lozenges. 

partition we can put into correspondence a q— weight g' 77 '. Notice that these q are not in 
any way related to the g-deformation parameter of Section II. The sum of g— weights of all 
plane partitions contained in a box r x s x t is known as q— enumeration of plane partitions 
or the partition function of Young diagrams 

z q{ r, M)=E? kl =nnii _ L 3 - + *- 2 = n n x • (32) 
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This formula is MacMahon generation function for the boxed plane partitions [2(| 
the q — > 1 limit this formula gives the number of the plane partitions containing in a box 
r x s x t. 

Let us consider the scalar product of the state vectors (j2HI) and (|3T)j) : 

N N 

S(N, M\{v}, {u}) = (0| n C( Vj ) II B( Uj )\0), (33) 

i=i j=i 

where {«} and {f } are the sets of independent parameters. We shall show that after the 
parametrization 

Uj = q 2 , Vj = q 2 , (34) 

the scalar product (}3*3*j) is the generating function (}3*2"|) of the plane partitions containing in 
a box N x N x M: 

J-l 



g— (0| n C(qi)Y[B(q^)\0} = Z q (N,N,M). (35) 
j'=-i i=i 

To make a connection of the scalar product (|33|) with the problem of enumeration of 
boxed plane partitions we shall use a graphic representation of its elements introduced in 
the previous Section. Consider a two-dimensional square lattice with 2N x (M + 1) sites. 
First N vertical rows of the lattice are associated with the operators C(vj) and the last TV 
vertical rows with operators B(uj) (see FIG. 1 and FIG. 2). The horizontal rows of the lattice 
are associated with the local Fock spaces, z-th raw with the i-th space respectively. Each 
horizontal edge of a lattice is labelled then with an occupation number rij of a correspondent 
Fock vector \rij)i. The scalar product (}3*3*j) is equal to the sum over all allowed configurations 
on a square lattice with the arrows on the first iV vertical rows pointing inwards, on the last 
N ones pointing outwards; on the right and on the left boundaries all occupation numbers 
are zeros. 

These configurations may be represented in terms of the lattice paths. The allowed 
configurations are then the number of possible non-crossing N paths starting from the down 
left N lattice edges (— iV;0),(— iV + 1; 0), (— 1; 0) and ending at the iV top right ones 
(1; M), (2; M), (N; M). The m-th path is running from (-N + m - 1; 0) to (to; M), 
1 < m < N. In the vertical direction paths follow the arrows, and only one path is allowed 
on a vertical lattice edge but any number of paths can share the horizontal ones. The number 
of paths sharing a horizontal edges is equal to the corresponding occupation number of the 
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edge. The length of the paths is (N + M). One of the possible configurations is represented 
in (FIG. 6). 

(-4,5) (-3,5) (-2,5) (-1,5) (1,5) (2,5) (3,5) (4,5) 
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r 



(-4,0) (-3,0) (-2,0) (-1,0) (1,0) (2,0) (3,0) (4,0) 
FIG. 6: A typical configuration of admissible lattice paths. 

The cells of the lattice under the m— th path may be considered as diagram of corre- 
sponding partition and may be thought of as the m— th column in the array (71^). The 
configuration of the paths in (FIG. 6) corresponds to the plane partition in (FIG. 4) and 
respectively to the array (p?T]) . 

On the other hand we can associate the vertical and horizontal edges carrying paths with 
lozenges. Lozenge (a) in (FIG. 5) corresponds to a vertical line of the path, while lozenge 
(c) to a horizontal one. Lattice edges without the paths correspond to a lozenge (d). The 
lozenge tiling is simply the projection of three-dimensional Young diagram with gradient 
lines. This establishes the mapping of the configurations generated by the scalar product 
(JHHJ) on the plane partitions. 

Consider the scalar product Due to commutation relations it is a symmetric 

function of iV variables Vj and also a symmetric function of N variables Uj. It is easy to 
verify that the number of operators C(u) should be equal to the number of operators B(u), 
otherwise the scalar product is equal to zero. The scalar product is evaluated by means of 
commutation relations (J20"j) - (J22j) . In the simplest case (N = 1) the scalar product is equal 
to 

S(l,M\v,u) = (0\C{v)B{u)\0) 

= g(v,u) {a(v)d(u) — a(u)d(v)} 
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u\ M+1 fu^- M - r 



x 



V J \V J J u _ I u 



-1 ) 



where g(u, v) is the element of i?-matrix (fT^|) and a(u) and d(u) are the eigenvalues of A{u) 
and D(u) (gj). 

For the arbitrary iV one may get Q, Q : 

2V 2V 

wMiw,^}) = (oinc(^)n 5 K-)i°) 

n^-^rnf^-^)' 1 }^ (36) 

where the matrix elements of N x N matrix H are equal to 

{/ \ M+N / v -M-N) . 

The parametrization (JMj) transforms the scalar product into 

S(N,M\{q}) = (-1)^ (n (g 2 ^ -g-^) _2 }detH, (38) 

where 

fc+j-i fc+j-i 
q 2 — q 2 



with s = q M+N . The determinant of the matrix Ti was considered in 
the alternating sign matrices enumeration problem and is equal to 

JV(JV-l) J T-r / t±. ^t±\ 2 I T-r S^q^ — 



221 in connection with 



detw=(-i)^<n(^-<r^j n ^ • ( 4 °) 

,i>& J i<j,k<N q 2 — q 2 



Therefore, 



S(N,M\{q}) = q-^~ J] 

1 „■ 1-^- AT -L 



1 _ gN+M+j-k 



jk+j-1 



l<j,k<N x 9 

= « 2 n • ( 4i ) 

l<j,k<N L y 

Finally, we have obtained the equality (J55J) for the scalar product: 



(0|C( ? -^)...C(g- 1 )C(g^) J B(l) J B(^)... J B(g^)|0) = g-^Z 9 (iV,iV,M). (42) 
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V. COORDINATE FORM OF STATE VECTORS AND SCHUR FUNCTIONS 



Using the explicit form of the operators B(X) we may rewrite the iV— particle state vector 
(129)1 in the "coordinate" form 



N 



M 



|*jv(mi,w 2 , -,u n )) = Y[ B(u k )\0) 



k=l 



Yl f{n}(u 1 ,u 2 ,... 1 u N )Y[\n j ) j , 

3=0 



(43) 



0<tiq ,71-^ , . . . ,n ^ <iV 
no+n 1 + ...+n M = N 



with the function / equal to 

f{n}(uuU 2l ...,U N ) = f( njli n j2t ... i n jM )(Ul,U 2 , - 9 U N ) = 



u 



N 



(44) 



Here the sum is taken over all admissible iV paths with paths starting from (l;ji), rij 2 
from (1; j 2 ), and n j&1 from (l;j M ) respectively, j l > j 2 > ... > j M ; n jl + ...+n jM = N, n jk ^ 0. 
The power tf is equal to the number of the d(u) vertices in the k— th vertical line of the 
grid, while t% is equal to the number of a(u) vertices respectively. 



a 



1 2 3 4 5 




FIG. 7: Lattice paths representation of a particular term of the 5-particle state vector, and a 
correspondent part of a plane partition. 



The following representation is valid 



/{„}(«!, u 2 , ...,u N ) = (uxu 2 ■ ... • u N ) S{\y{ul,v%, ...,u 2 N ) 



where S{\}(ui,u 2 , ...,un) is the Schur function 



S(^x 1 ,...,x N )(ui,u 2 , ...,u N ) 



det «" i+Al ) 



(45) 



(46) 



Y\l<i<j<N\ u i ~~ u j) 

and there is one to one correspondence between the occupation number configuration {n} 
and partitions {A}: 
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{n} = (n jl ,n j2> ...,n jM ) {A} = (Ai,..., X N ) 
Afc = ii; 1 < fc < 
Afc = J2; n h + 1 < fc < n j2 + n^j 



(47) 



Afc = i/v; + ... + n h[ _ 1 + 1 < k < n h + ... + n jM = N. 



The state conjugated to (|4*3^) is 



N 



M 



(^ N (u u u 2 , ...,u N )\ = (0| [] C(u k ) 

k=l 



E /wK 1 ,^ 2 ,-,^ 2 ) II( n ili' ( 48 ) 



0<nQ ,. .. ,n^y- < A/" 
™0+ ri l + --- + ri M= ^ 



i=o 



with 



t?-t? d-t a 



"l U 2 ■•• "TV 1 



(49) 



where the sum is taken over all admissible N paths with nj 1 paths ending at (— l;ji), rij 2 
ending at (-1; j 2 ), and n jM ending at (1; j M ) respectively, j t > j 2 > ... > j M ] n h + ... + n jM = 
N, n jk + 0. 



-5 -4 -3 -2 -1 



FIG. 8: Lattice paths representation of a particular term of the conjugated state vector. 



Due to the orthogonality of the Fock states (fTTj) the scalar product is equal to 
S(N,M\{v},{u}) = (^ N (v u v 2 , ...,v N )\^ N ( Ul ,u 2 , ...,u N )) 

= J2 f{n}(vT 1 ,V2 1 ,-,V^ 1 )f{n}(ui,U 2 ,...,U N ) 



0<riQ ,...,n^j<N 
nQ+n 1 + ...+n M = N 

j.d mcl j-d j.a 



0<n ,n 1 ,...,n M <N C 
^0+ n l+---+ n M~ N 



B 



j-d _j-a 

u N 



(50) 
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Taking into account the representation ()45|) we can express the scalar product in terms of 
Schur functions 



/ N M \ 

S{N,M\{v},{u})=[U-* s ) E ^A}^,^,...,^)^}^ 2 ,^ 2 ,...,^ 2 ), (51) 

V=l j I AC{Af^} 

where the sum is over all partitions, A , into at most N parts each of which is less than 
or equal to M. Comparing this formula with (|36|) we obtain the following determinantal 
expression 

E S {X }(xl,xl,...,x 2 N )S {x} (yl,y2,... J y 2 N ) (52) 

XC{M N } 

\j=l ) \ j>k Uj t>k l> m x l x m ) 

where the matrix H is (137(1 

H jk = {(xkVj) M+N ~ (x k yjy M ~ N } x 



XkVj - {xkVj) 1 

The volume \tt\ of the plane partition ir in a box (N x N x M) may be expressed as 

-N N 

2|tt| =N 2 M+Y: k (Zj - It) + EO' - 1) 0; " , (53) 
fc=-i i=i 

where the first sum in this equality is over the columns going along the " s" side of hexagon 
while the second one is over the columns along the "r" side respectively, and l^ d is the 
number of lozenge of type a, d (see FIG. 5) in the j-th column of the hexagon. It may be 
checked that lj — lj = tj—tj. The substitution of the parametrization (|33|) into the equation 
(IHOJ) gives for the scalar product 

S(N,M\{q}) = (^ iV (g-^ ? - 1 ,...,g-f)|vl> iV (l, g ,...,gV)) 
= q ^ E? ■ 



Together with (J4*T|) this equation provides us with the proof of MacMahons enumeration 
formula for the boxed plane partitions within the frames of Quantum Inverse Method. 
From the equation ((5T|) we obtain the equality 



E ?, g*- 1 )^, g 2 , ...,^)=E ? W - 

AC{A/ JV } P-P 
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The parts of the partition A (jTfj) may be considered as the coordinates of the particles 
(the coordinate A^ corresponds to the k-th particle). The state of the system is spanned by 
the orthonormal basis |A) = |Ai, A2, A m ) : (/x| A) = 5{^},{a}- The iV-particle state vector 
(jIHJ) is given then by the equation 

\^n(u 1 ,u 2 , -,u n )) = (miu 2 ■ ... ■ u N y M s {X}(ul,ul, ...,u%)\\), (54) 

\C{M N } 

and respectively 

JV 

Y[B( Uj )\0) = J2 S {x} (ut,u 2 2 ,...,u%)W, (55) 

3=1 XC{M N } 

where the sum in both formulas is taken over all partitions, A , into at most N parts each 
of which is less than or equal to M. By the construction operators B{u) = u~ M B(u), and 
C{u) = u M C{u) possess the following properties 

B(l)|A> = ElA*>, <A|C(1) = 53<Ml, (56) 

/OA fiDX 

and 

N 

Y[B( Uj )\X) = ^S W x } (uI,i4,---,u 2 n)\v), (57) 

3=1 /OA 
N 

(AIJJC'K) = Yl S MX }(uf,U2 2 ,...,u^ 2 )(fi\, 

3=1 /OA 

where S{n/\}(xi, x^, xn) is a skew Schur function indexed by a pair of partitions fi and A 
such that A C fi |21|. From this point of view operators B{u) and C{u) may be considered 
as the transition operators between the diagonals of a plane partitions. 



VI. CORRELATION FUNCTIONS 



Let us calculate the generating function of the plain partitions contained in a box N x 
N x M provided that the height 7Tin of the stack of the cubes is fixed and equal to m. To 
find this function we have to consider a scalar product on a lattice under the condition that 
the iV-th lattice path enters the N-th vertical line of the grid at the m-th row: 

N N-l 

P 1N (m\{v}, {u}) = (0| I] C( Vj ) I] B( Uj )<f>l\0). (58) 

3=1 3=1 
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1Z N ] 




FIG. 9: A Young diagram with fixed heights it of the stack of cubes. 

To calculate this scalar product we may use the following decomposition of the operators 
B{u) and C(u): 

M 



B(u)\0) = u- M J2<py>\0), 

j=0 

M 

.M/n\ \ ^ a „.-2j 



(0\C(u) 



U 



uu 



(59) 
(60) 



3=0 



The substitution of decomposition ()59|1 into the scalar product gives 



M 



S(N,M\{v},{u})=J2( u n 



-M+2m 



P 1N (m\{v},{u}). 



(61) 



m=0 



The determinant of the matrix H in (J36j) may be developed by the last column. The 
comparison of the obtained decomposition with (|5T|l leads to the equality 
P 1N {m\{v},{u}) 



f n MM^t) 2 ) detQ 

S(N,M\{v},{u}) M^lli (u t ) 2 



detH' 



The entries of N x N matrix Q are given by 



QjN = (Vj) 



M+N-l-2m 



, Qjk = Hjk, k N, 



where Hjk are the matrix elements (JHTj) . After the parametrization 
probability of the height tvin to be equal to m is 

PiN(m\{q}) 



(m) 



IN 



S(N,M\{q}) 



(62) 
(63) 

we find that the 
(64) 



n(i-9 



N-t-2 



*t<N 



detQ 
detH' 
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where 



j(M+N-l-2m) 

QjN = q a , Q ifc = H jk , k^N, (65) 



and Tijk are the matrix elements ()H9j) . It is evident that this expectation value is the same 
for the height uni in the opposite corner of the diagram. 

The correlation function of the heights of the columns i^ini^ni at the opposite sides of 
the Young diagram (see FIG. 9) may be obtained from the following scalar product 

N-l 2V-1 

PNi;iN(n;m\{v},{u}) = (0\<j> n I] C( Vj ) J[ B( Uj )<i>l\0). (66) 

j=i j=i 

The scalar products (JoTj) and (JSfij) may be expressed in terms of the skew Schur functions 
(JSZJ) as well. 

The other function of interest is the projection of the Bethe wave function ITjli^KOlO) 
on the "steady state" vector: 

M 

\P)= E IlK-)r (67) 

0<tiq ,ni , .. -,riM <N j=Q 
n 0+ n l+---+ n M =N 

From the representation ()43|) and relation we obtain the equality 

N 

(P\HB{ Uj )\0) = { Ul u 2 -...-u N y M Yl S {x} (ul,ul...,u 2 N ), (68) 

i=l AC{M W } 

where the sum is over all partitions, A , into at most N parts each of which is less than or 
equal to M. It is known that j^j] 



det (x{ 1 



M+2N-J 



E S { x}(x 1 ,x 2 ,...,x N ) = — — T — 2jV _ A (69) 



n r 



i=l ^ l<i<j<N 

and we obtain for the projection ()68jl 

AT AT 2 ^ 



(^in^)io> = n T ^ ? n r _^ 2 - 

i=l i=l L u i l<i<j<N 1 a i a j 



(70) 



By the construction the considered correlation function is the generating function of the 
symmetric plane partitions, the plane partitions satisfying the condition Hij = n^. From 
flSJ) and the relation tf - t% = l{ - 1% it follows that 

N 

(pin s («f)i°>=E«? 



id ja id la id la 

,1 In. 2 l 2 „,N l N 

a 2 ' ■■■ ' U N 



i=i 
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The volume of the symmetric plane partitions may be expressed as 



N 

2\ir\ syrn = N 2 M + £(2j - 1) (zj - if , (71) 
where the sum is over the columns going along the "r" side of the hexagon. Then 

N 2 
j=l sym 

and we obtain the well known result for the generating function of the symmetric plane 
partitions 

E 5 {A} (g,g 3 ,...,^- 1 ) = E^ 1 - 

\C{M N } sym 

Till now we have considered plane partitions in a box with r — s. To study the general 
case when r^swe have to consider the following scalar products 

N L N 

S A (N,L,M\{v},{u},{u A }) = (0|nCK)n^«)n 5 (%)|0), (72) 

j=l j=i j=i 

N L N 

S D (N,L,M\{v},{u},{u D }) = (0|nC(^)II^«)n^)|0)- (73) 

j=l j=\ j=l 

Following the mapping introduced in this Section it may be shown that these averages are 
the generating functions of plane partitions in a box r x s xt with r = N,s = N + L,t = M 
for ((ZD) , and r = N + L, s = N, t = M for (Zfl . 

VII. BOXED PLANE PARTITIONS AND TODA LATTICE 

The N x N matrix (J39|) is a Hankel matrix with the matrix elements 

H jk = h{ q 3 +k - v ), (74) 

where 

M+N M+N 

h(q) = q \ ~ q \ 2 . (75) 
qz — q 2 

By the successive subtraction of columns the determinant of this matrix may be brought 
into the form 

detH = deth, (76) 
h jk = Dl +k - 2 h(q), 
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and D q is the g-difference operator: 

D q f(z) = f(qz)-f(z). (77) 



Following the standard procedure 
det Ti. satisfies the equation 



251 ] it may be shown that the function t(N, M; q) 



which, after the substitution 

.„ . , T(N + l,M;q) 
"W"' = 1 " r { N,M- q) 

becomes the g-difference Toda equation: 

D 2 p(N; q) = e P(N+l;q)-p(N;q) _ e p(N;q)-p(N-l;q) _ 

The role of time plays the deformation parameter q. 
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